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Abstract. We present an evolutionary sequence of a low mass 
star from the Asymptotic Giant Branch (AGB) through its post- 
AGB stage, during which its surface chemical composition 
changes from hydrogen-rich to strongly hydrogen-deficient as 
consequence of a very late thermal pulse, following the so- 
called born-again scenario. The internal structure and abun- 
dance changes during this pulse are computed with a numeri- 
cal method which allows the physically consistent calculation 
of stellar layers where thermonuclear and mixing time scale are 
comparable — a situation which occurs when the helium flash 
driven convection zone extends to the hydrogen-rich surface 
layers during the pulse peak. The final surface mass fractions 
are [He/C/O]=[0.38/0.36/0.22], where the high oxygen abun- 
dance is due to diffusive overshoot employed during the AGB 
evolution. These models are the first to achieve general agree- 
ment with the surface abundance pattern observed in hydrogen- 
deficient post-AGB stars — e.g. the PG 1 159 stars or the WR- 
type central stars of planetary nebulae — , confirming the born- 
again scenario with a physically consistent calculation and sup- 
porting the occurrence of convective overshooting in thermally 
pulsing AGB stars. 
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1. Introduction 

Stars on the so called Asymptotic Giant Branch (AGB) have 
strong stellar winds, which gradually reduce the mass of the 
hydrogen-rich stellar envelope. When this envelope mass falls 
below a critical value, the stars leave the AGB to become post- 
AGB stars, central stars of planetary nebulae (CSPNe), and fi- 
nally white dwarfs. Post-AGB stars show a variety of surface 



abundances ( |Mendez, 199T| ). About 80% of all CSPNe show a 
solar-like composition while the remaining ones are hydrogen- 
deficient. Among the latter are Wolf-Rayet type CSPNe ([WR]- 
CSPNe) and the extremly hot PG 1159 stars with typical sur- 
face ab undan ces of [He/C/O]=[0.33/0.50/0.17] (Dre izler and 
Heber (1998); see also Koesterke and Hamann (1997) and ref- 



erences in both papers). Hydroge n-deficiency is also found i n 
white dwarfs of spectral type DO ( Dreizler and Werner, 1996 ). 



The origin of the hydrogen-deficiency in post-AGB stars 
is a longstanding problem. Most post-AGB calculations pre- 



diet a hydrogen-rich surface composition ( 


Schonberner, 1979 


[Schonberner, 1983; 


Wood and Faulkner, 1986; Vassiliadis and 


Wood, 1994; 


Blocker, 1995a; 


Blocker and Schonberner, 1997). 



So far, no post-AGB models reproduced the observed high 
carbon and oxygen abundance. The most promising scenario 
for obtaining a hydrogen-defic ient surface com position en- 
vokes a very late thermal pulse ( Fujimoto, 1977 ; Schonberner, 
1979; Iben et al., 1983 ) — i.e. a pulse which occurs af- 
ter the star has already left the AGB — during which the 
pulse driven convection zone can mix hydrogen-free mate- 
rial out to the stellar sur face. Within this born-again sce- 
nario, Iben and McDonald (1995) obtain surface mass fractions 



of [He/C/O]=[0.76/0.15/0.01], i.e. their model indeed became 
strongly hydrogen-deficient. Hovever, the large oxygen abun- 
dance found in most H-deficient post-AGB stars could not be 
reproduced be these, nor by any other calculation. These diffi- 
culties have posed a strong limitation to the whole scenario. 

In this Letter we present a post-AGB model sequence start- 
ing from an AGB model computed with overshoot (Herwig 
et al., 1997), and using a numerical method of computing nu- 
clear burning and time-dependent convective mixing simulta- 
neously. 



2. Numerical method 

The stellar models are based on the stellar evolution code 
described by Blocker (1995b). However, the treatment of 
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the chemical evolution was entirely replaced by a numerical 
scheme which solves the time dependence of the considered 
nuclear species — i.e., the changes due to thermonuclear reac- 
tions and due to mixing — in one single step. This enables us, 
in contrast to earlier investigations of very late thermal pulses, 
to reliably predict the chemical abundance profiles and the nu- 
clear energy generation rates in situations where the time scales 
of nuclear burning and mixing are comparable. The abundance 
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change for each isotope at each mesh point due to diffusive 
mixing and nuclear processing is given by 





d 




dm 



om 



(1) 



where Xj contains the abundances of all considered isoptopes 
at the j*'' mesh point, Fj is the nuclear rate matrix, D is the dif- 
fusion coefficient describing the efficiency of convective mix- 
ing, r is the radius, m the mass coordinate and p the density. 
This leads to a set of non-linear equations with M ■ N un- 
knowns, where M is the number of grid points and N is the 
number of isotopes. In the present calculations, M is of the or- 
der of 2000, and = 15 as the main thermonuclear reactions 
for hydrogen burning through the pp chains and the CNO cy- 
cle as well as the main helium burning reactions are included. 
The solution is obtained fully implicit with a Newton-Raphson 
iteration scheme by making use of the band-diagonal structure 
of the problem. The scheme converges to sufficient precision 
within about 3 iterations. A coupled solution of one nuclear 
reaction at a time and time-dependent mixing, including also 
the structure equations, has already been applied by Eggleton 



(1972) 



O 



3. The AGB starting model 

We start with an AGB model with A/zams = 2 Mq which 
has been evolved over 22 thermal pulses, including convective 
overshoot at all convective boundaries. The treatment and ef- 
ficiency (f=0. 016) of overshoot is the same as in Herwig et 
al. ( 1997 ). In comparison to models without overshoot the in- 
tershell region is much stronger enriched in carbon and oxygen 
(mass fractions [He/C/O ]=[0.35/0.43/0. 19]) , which causes a 
stronger third dredge-up ( |Herwig et al., 1999| ). At the 16*'' ther- 
mal pulse (TP) the hydrogen-free core has a mass of Afcorc ~ 
0.573 M0 and dredge-up starts to operate, leading to a carbon 
star model at the last computed TP. At this stage the model star 
has a total mass of M = 1.42 Mq and Mcorc = 0.604 Mq. We 
then artificially increase the mass loss (M > 10^^ M©/ yr) in 
order to force the model to leave the AGB at the right phase 
to develop a very late thermal pulse. This procedure is justified 
for this exploratory work because it does not affect the nucle- 
osynthesis and mixing during the very late TP. 

4. Evolution through the very late thermal pulse 

Two cases of the born-again scenario should be distinguished. 
Depending on the time when the post-AGB thermal pulse oc- 
curs, shell hydrogen burning may still be active or may already 
have ceased. 

In the first case, the He-flash driven convection zone can- 
not extend into the hydrogen-rich envelope due to the entropy 
barrier generated by the burning shell (Iben, 1970- In the sec- 
ond case, which is realized in our model sequence, hydrogen 
shell burning is extinct and the star has entered the white dwarf 
cooling domain (Fig.jl]). We designate a TP in this situation as 
a very late TP. 




Fig. 1. Track in the HR diagram of our post-AGB model dur- 
ing the evolution through a very late TP. At the first mark along 
the track, the He-flash has already caused a prominent convec- 
tively unstable region. At the second mark, the outwards grow- 
ing convection zone has reached the envelope and protons start 
to enter the convective zone (compare Fig. g and Fig. At the 
third mark the hydrogen luminosity has reached its peak (see 
Fig.^ and Fig.^. The surface composition is hydrogen-free 
beyond the last mark. 
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Fig. 2. Evolution of the top boundary of the convection zone 
(dash-dotted line with filled circle for every second stellar 
model, left scale) as it extends into the hydrogen-rich envelope, 
and the nuclear luminosity due to hydrogen burning (solid line, 
right scale). The four grey dots along the solid line correspond 
to the marks in Fig. |l]. The dashed line shows the mass coordi- 
nate of the hydrogen-free core which is for i > yr identical 
with the total stellar mass. 



As the helium luminosity increases in the course of the 
He-flash in our model sequence (first mark in Fig. |]), the cor- 
responding region of convective instability enlarges (Fig.^. 
When the upper convective boundary reaches the mass coordi- 
nate where the hydrogen abundance increases, convective mix- 
ing transports protons downwards into the hot interior (Fig. |]). 
The protons are at some point captured by via the reac- 
tion ^^C(p, 7)'^'^N. The peak of the resulting luminosity due to 
hydrogen burning (see also Fig.0) is located at the mass co- 
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Fig. 3. Internal structure and composition at the onset of hydro- 
gen ingestion into the He-flash convection zone during the very 
late TP (second mark in Fig. |]). Top panel: The nuclear energy 
generation is dominated by processing of helium at the bottom 
of the He-flash convection zone. Energy due to proton capture 
is released in the upper part of the He-flash convection zone 
(left scale). The hydrogen profile reflects the simultaneous nu- 
clear burning and convective mixing (right scale). The surface 
composition at this stage is still hydrogen-rich. Bottom panel: 
The diffusion coefficient (left scale) visualizes the convectively 
unstable region corresponding to the He-flash, which has just 
reached the lower part of the hydrogen-rich envelope. 



ordinate where the nuclear time scale equals the mixing time 
scale (~ one hour). 

The profile of hydrogen in Fig.|^ and ^ demonstrates that a 
coiTect treatment of simultaneous burning and convective mix- 
ing is essential for this evolutionary phase. A treatment of con- 
vective mixing which does not include the simultaneous com- 
putation of the isotopic abundances according to the equations 
of the nuclear network would fail to predict a correct hydro- 
gen profile. In particular, such a treatment would possibly let 
the protons travel too deep into the convective region, without 
considering that they would have been captured already on the 
way. Then, the energy generation rate due to proton captures 
may be overestimated and not correctly located. 

The energy from proton captures is released in the upper 
part of the He-flash driven convection zone, which leads to a 
split (at TOr = 0.595 Mq) of the convective region (Fig. The 
two convective regions are then connected by the overlapping 
overshoot extensions, but Fig. ^ shows that the second convec- 
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Fig. 4. Internal structure and composition at the time of max- 
imum energy generation due to proton captures (third mark in 
Fig.[l|). Top panel: Nuclear energy generation due to hydro- 
gen burning and helium burning (left scale) and hydrogen pro- 
file (right scale). Bottom panel: The diffusion coefficient (left 
scale) shows that the convectively unstable region of the ther- 
mal pulse is split into two (toi- = 0.595 M0). 



tive zone is only short lived since the amount of hydrogen avail- 
able in the envelope is quickly consumed. 

Figure Q shows that the hydrogen burning convection zone 
extends over ~ 10^^ M© and reaches from m,- — 0.595 Mq 
up to the surface of the stellar model. The surface hydrogen 
abundance declines rapidly due to mixing and proton captures 
in the deeper layers. The period of the largest hydrogen burning 
luminosity (shown in Fig.0) of in — 10* L© lasts for less 
than a week, and the whole episode of convective hydrogen 
burning is a matter of about a month. Overall, 5 • lO"'"^ Mq of 
hydrogen are burnt. At peak hydrogen luminosity the hydrogen 
mass fraction at the surface is 3.4 • 10^^ and the total amount of 
hydrogen still present in the star is Mj^ ~ 7.8 • 10^^ Mq. Thus, 
in this sequence the star is already hydrogen-deficient before it 
returns to the AGB domain in the HRD. 

Figure ^ shows abundance profiles before and after the 
mixing and burning event due to the very late TP. While 
the star still shows the typical hydrogen-rich AGB abun- 
dance pattern before the convective region has reached into 
the envelope (top panel. Fig. |]), the mixing during the con- 
vective hydrogen burning leads to a hydrogen-free surface 
with [He/C/O]=[0.38/0.36/0.22] and a mass fraction of 3.5% 
of neon. The step in the abundances of ^He, ^^C and ^^O at 
0.596 Mq (lower panel) coiTesponds to the split of the convec- 
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Fig. 5. Chemical profiles (mass fraction vs. mass coordinate) 
of the upper mass region also covered by Fig.^ Top panel: 
Profile before the convective region of the He-flash has reached 
the envelope corresponding to position of first mark in Fig. |l] 
Bottom panel: Profile corresponding to position of last mark 
in Fig.|l] after all hydrogen has been processed (only trace 
amounts of 4 ■ 10^^^ Mq are left in the whole star), the abun- 
dance at the surface differs only slightly from the intershell 
abundance. 

five region due to hydrogen burning. While the hydrogen burn- 
ing leads not to a significant abundance changes for the major 
isotopes (only 5 ■ 10~^ Mq of hydrogen are processed), helium 
burning continues to process helium at the bottom of the He- 
flash convective zone. The final surface abundances are very 
similar to the intershell abundances during the thermal pulse. 

After most of the hydrogen is burnt, the corresponding 
upper convection zone disappears when the local luminos- 
ity drops. It takes about one year until the He-flash convec- 
tion zone has recovered to its original extent (Fig.^. The star 
then follows the evolution as known from the born-again sce- 
nario (for a recent account on this scenario see Blocker and 
Schonberner, 1997). Energetically, the return into the AGB do- 
main is almost exclusively driven by the energy release due to 
helium burning, which exceeds the additional supply of energy 
from hydrogen burning by orders of magnitude. 

5. Conclusions 

Using a numerical method to treat nuclear burning and mix- 
ing simultaneously in stellar evolution calculations, which al- 
lows a reliable and robust modelling of very late thermal pulses, 
we have shown that the general surface abundance pattern ob- 
served in hydrogen-deficient post-AGB stars can be explained 
within the born-again scenario. Our new post-AGB sequence 



shows that due to the energy generation and convective mixing 
during a very late thermal pulse a born-again star forms which 
displays its previous intershell abundance at the surface. 

We have based the calculation on an AGB model sequence 
computed with overshoot, which shows a high carbon and oxy- 
gen intershell abundance. Thus, the fact that the abundance 
pattern of our post-AGB model after the thermal pulse agrees 
with the observation of hydrogen-deficient post-AGB stars like 
PG 1159 and [WC]-CSPNe strongly supports the assumption 
of extra mixing beyond the convective boundary of the He- 
convection zone in AGB stars. We conclude that the very late 
thermal pulses can indeed be identified as one cause for the 
hydrogen-deficiency in post-AGB stars. 

However, we note that not all H-deficient post-AGB stars 
are completely free of hydrogen (Leuenhagen and Hamann, 
1998), as predicted by our model. Other possibilities than the 
born-again scenario to achieve H-d eficiency might also exist 
( fTylenda, 199^ ; [Waters et al., 1998[ ). Whether post-AGB mod- 
els which are not entirely hydrogen-free can be obtained within 
this scenario requires a study of the variation of the late thermal 
puis with the inter-pulse phase at which the star leaves the AGB 
(Iben, 1984), and possibly the consideration of other mixing 
processes, e.g. due to rotational effects (Langer et al., 1999), 
which has to be left to future investigations. 
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